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Разработана и реализована математическая модель нестационарных гидравлических процессов в промышленных центробеж- 
ных каскадах обогащения урана, предназначенная для применения в качестве экспертной системы в автоматизированной систе- 
ме управления технологической схемой и использования в компьютерном тренажере для подготовки специалистов раздели- 
тельного производства. 


Введение 

Автоматизированная система управления тех- 
нологической схемой (АСУТС) центрифужного 
разделительного завода предназначена для обеспе- 
чения и поддержания параметров расчетной техно- 
логической схемы (ТС) при изменяющихся внеш- 


них условиях. АСУТС основана на моделировании 
стационарных процессов разделения. Многолет- 
ний опыт эксплуатации показал её эффективность, 
особенно в последние годы, когда существенно вы- 
росли требования по качеству выпускаемой про- 
дукции. 
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Вместе с тем, большое влияние на работу разде- 
лительного производства оказывают нестационар- 
ные гидравлические процессы, возникающие как в 
ходе нормальной эксплуатации каскадов, например, 
при изменении частоты питающего тока, при пере- 
ходе на новый технологический режим, при плано- 
вом отключении части оборудования и т.п., так и 
при нарушениях нормальной работы схемы, напри- 
мер, при закрытии внешних потоков каскадов. 

Эти процессы влияют на качество выпускаемого 
товарного продукта, эффективность работы ТС за- 
вода и могут привести к недопустимым перегрузкам 
оборудования. Во всех нестационарных процессах 
необходимо обеспечить безопасность оборудования 
и свести к минимуму потери работы разделения, что 
и обуславливает важность изучения и учета этих 
процессов. Кроме того, по сравнению с диффузион- 
ными, центрифужные каскады имеют меньшее га- 
зосодержание, а следовательно и меньшую инер- 
ционность. Малая инерционность приводит к уси- 
лению влияния нестационарных процессов на эф- 
фективность работы разделительных каскадов. 

В связи с этим возникла необходимость реше- 
ния задачи полномасштабного учета нестационар- 
ных гидравлических процессов в разделительном 
производстве. 

С целью решения поставленной задачи разрабо- 
тана и программно реализована математическая мо- 
дель нестационарных гидравлических процессов в 
многокаскадной ТС разделительного производства. 

В отличие от ранее известных работ [1], рассма- 
триваемая модель нестационарных процессов учи- 
тывает особенности построения ТС разделительных 
предприятий, в частности, заключающиеся в том, 
что соединения ступеней в каскаде могут осущест- 
вляться по параллельно-последовательной схеме. 

Кроме того, в модель включены линии межка- 
скадных коммуникаций (МКК), связывающие ка- 
скады между собой. 

1. Описание математической модели 

нестационарных гидравлических процессов 

1.1. Объемы, входящие в модель 

В описываемой модели газосодержание ступе- 
ни п каскада т многокаскадной ТС считается со- 
средоточенным в 10-ти объемах (для 8-секционной 
ступени). Из них 8 объемов роторов газовых цен- 
трифуг (ГЦ) секций і с газосодержанием М 1іт сту- 
пени п каскада т (і,п,т - номер секции, ступени, 
каскада). Газосодержание М ІІШ определяется давле- 
нием на носике отборника тяжелой фракции после 
прохождения скачка уплотнения Р Мш [2]: 

М\іпт — Ѵц т т (® )-^ 0 іпт > — 1 — 8 , 

где М ІІШ - газосодержание роторов ГЦ /-ой секции 
л-ой ступени каскада т; Ѵ ІІШ - коэффициент, зави- 
сящий от частоты вращения роторов со и типа ГЦ. 


Аналитическое выражение для газосодержания 
М Шт определено исходя из принятой упрощенной 
газодинамической модели ГЦ. 

В газовые центрифуги секции і ступени п пода- 
ется поток питания 1 Шт а выходят потоки легкой 
Ьлітт и тяжелой фракции Ь Піт . 9-ый объем трассы 
отвала ступени п каскада т с газосодержанием М 1ш 
и давлением Р Тт включает в себя коммуникации от 
коллекторов отвала отдельных ГЦ до межступен- 
ного регулятора давления. Его газосодержание 
определяется выражением: 

М = V Р 

где М 1т - газосодержание трассы отвала ступени п 
каскада м; Ѵ 1ш — константа, зависящая от физиче- 
ских свойств рабочего газа и геометрических раз- 
меров трассы отвала ступени п каскада т; Р Тт - да- 
вление в коллекторе отвала ступени п каскада т. 

В данный объем входят секционные потоки 
отвала Ь Піт , а выходит суммарный поток отвала 
ступени С Тт равный расходу О рт через межступен- 
ный регулятор. Принято, что поток в секционной 
(секция І) трассе отвала Ь ПІШ равен потоку отвала из 
ГЦ секции і Ь піш . 

10-ый объем трассы питания ступени п каскада 
т с газосодержанием М іш и давлением Р 0ш может 
включать следующие коммуникации: 

• блочная трасса отбора предыдущей ступени; 

• байпасная линия отбора предыдущей ступени; 

• коллектор питания ступени; 

• байпасная линия отвала следующей ступени; 

• блочная трасса отвала следующей ступени. 

Для газосодержания трассы питания ступени п 

каскада т имеем: 

М =Ѵ Р 

х Ъпт г Ъпт л Опт 1 * * * 5 

где М 3т - газосодержание трассы питания ступени 
п каскада т; Ѵ 3т - константа, зависящая от физи- 
ческих свойств рабочего газа и геометрических раз- 
меров трассы питания ступени п; Р І)ш - давление в 
коллекторе питания ступени п каскада т. 

В общем случае в этот объем могут входить сле- 
дующие потоки: 

• отбора двух предыдущих ступеней О л „_ 1т С Лп _ 2м , 

• отвала следующих двух ступеней 0 Тп+1т , С Тп+2 ,т- 

• питания ступени Т ш и Т Пю являющиеся пото- 
ками отбора и отвала полки к (кфт) соответ- 
ственно; 

• внешнего питания данной ступени Г 0 „Ц), яв- 
ляющегося заданной функцией от времени ?; 

• воздуха С в „(1) являющегося заданной функцией 
от времени Г при моделировании нарушения ва- 
куумной плотности оборудования. 

Из 10-го объема выходит сумма потоков пита- 
ния секций данной ступени Ь Шіш . 
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Технические науки 


1.2. Основные уравнения нестационарной гидравлики 


Основные уравнения нестационарной гидра- 
влики - это уравнения баланса вещества в выде- 
ленных объемах, которые не учитывают коррози- 
онные потери. Изменение газосодержания в каж- 
дом выделенном объеме определяется разностью 
входящих и выходящих потоков. Для выделенных 
объемов ступени п каскада т можно записать сле- 
дующие уравнения нестационарной гидравлики: 




А)1 іпт Аиш А'У Ііпт ’ * 1,2,. ..,8, (1) 


дМі пт _ . г рч 

й* /Ц, 7ТШет 17 рпт ’ 'А/ 


Й ІМ, 


ж 




&лп-2^Лп-2,т + ^Тп+ 1^*Тп+ 1,ет ~*~ 


"^^гп+г^ги+г.т + ^ор (^шт Т ОРт (/ ) ) 

“ХА» ш », (3) 

1 = 1 


где К- количество каскадов ТС; и - приз- 
наки связи между потоками отбора и отвала каска- 
да ^ с питанием ступени л (кф>п); 8 Лп _ ь 8 Лп _ ъ 8 ШЬ 
8 Тп + 2 , (5 Й , - признаки наличия потока (0 или 1); 
<5ол - признак отборной ступени (0 или 1); С Шт , 
ТоРт - соответственно потоки отбора отборной сту- 
пени и каскада т. 

Ур. (1) составляется для каждой секции ступени 
п каскада т. Уравнения (2) и (3) составляются для 
каждой ступени п каскада т. 


1.3. Уравнения для потоков 


Поток питания отдельной ГЦ Ц ит проходит че- 
рез расходомерную диафрагму. Для потока питания 
ГЦ секции / ступени п каскада т имеем: 


т = N К Р 

■НО 1 іпт 1 ѵ іпт Л 0 іпт Г 0іпт ’ 


(4) 


где М іш - количество ГЦ в секции і ступени л; К Шш - 
коэффициент расхода; Р Шт - давление на входе в ГЦ. 


Секционный поток питания Ь Шт для принятого 
турбулентного течения газа по трассе питания [3] 
определяется следующим выражением: 


(А* іпт ^>0 іпт 


= Р 2 - Р 2 

0 пт Ош/и 


(5) 


где Р 0ш - давление в коллекторе питания ступени; 
Р Шт - давление на входе в ГЦ; % 0іт - коэффициент 
гидравлического сопротивления [3]. 

Исходя из ламинарного характера течения газа 
в отборнике тяжелой фракции ГЦ, для потока 
отвала ГЦ секции / ступени л каскада т Ь Піт имеем: 


ТХіптГЭ Тіпт р 2 р 2 

~ Г \ѵ0 іпт Г Тпт ’ 


N.. 


( 6 ) 


Поток отбора ГЦ /-ой секции Ь ЛІІІІ определяем с 
помощью формализованной универсальной гидра- 


влической характеристики ГЦ, которая предста- 
вляется в виде: 



где Ь I Іт —Ъ6 тт - коэффициенты гидравлической ха- 
рактеристики, зависящие от типа ГЦ, температуры 
помещения, температуры охлаждения ГЦ, концен- 
трации легких примесей, частоты вращения ® и т.п. 

К уравнениям (4-7) в принятой модели добав- 
ляются уравнения, связанные с работой межсту- 
пенного регулятора. 

Считая истечение газа через сопло регулятора 
критическим, запишем расход через регулятор в виде 

Ср=0(И)Р до , (8) 

где В - коэффициент расхода, зависящий от поло- 
жения клапана регулятора А; Р до - давление до ре- 
гулятора. 

В данной модели рассматриваются пропорцио- 
нальные регуляторы с регулированием по перепаду 
давления Р Тп -Р 0 „. В этом случае линеаризованное 
уравнение регулятора можно записать через прира- 
щение давлений А Р Тп , АР {] „: 

+ Р = Пп(АР Тп - АР 0п ) , (9) 

где ц=ДА/А - относительное изменение положе- 
ния клапана регулятора; Тп и Пп постоянная вре- 
мени регулятора и коэффициент усиления регуля- 
тора ступени и; А - величина, характеризующая по- 
ложение регулирующего клапана регулятора. 

Внутренние потоки питания Т Мю Т Шп определя- 
ются через соответствующие потоки отбора и отва- 
ла каскадов с учетом транспортной задержки т, м , 
под которой понимается время прохождения ги- 
дравлического возмущения по соответствующим 
линиям МКК. Время т. аі рассчитывается исходя из 
уравнения неразрывности для установившегося од- 
номерного потока, изолированного от внешней 
среды [2], а также из расчетных величин потоков и 
давлений на участках линий МКК и заданных зна- 
чений геометрических размеров этих участков. 

2. Схема решения уравнений 

нестационарной гидравлики 

Расчет нестационарных гидравлических про- 
цессов сводится к нахождению величин давлений 
во всех введенных объемах для всех ступеней ка- 
скадов, положений регулирующего клапана меж- 
ступенных регуляторов, как временных функций 
при заданных начальных условиях. В качестве на- 
чальных условий приняты расчетные значения да- 
влений исходного стационарного состояния ступе- 
ней каскадов. 
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Входящие в уравнения (4-9) коэффициенты 
расхода шайб, коэффициенты гидравлических со- 
противлений трасс, коэффициенты гидравлических 
характеристик ГЦ, параметры регуляторов и т.п. 
предполагаются известными из экспериментальных 
и расчетных данных. Подставляя выражения для га- 
зосодержания и потоков в уравнения (1-3) получим 
систему дифференциальных уравнений первого по- 
рядка относительно давлений Р МІШ , Р Тт , Р ш , Р Шт и 
положений регулирующего клапана регуляторов ц т 
для каждого из каскадов ТС. 

Таким образом, задача моделирования неста- 
ционарных процессов в промышленных центри- 
фужных каскадах сводится к решению системы 
дифференциальных уравнений первого порядка с 
заданными начальными условиями. Задача Коши 
за редким исключением не имеет аналитических 
методов решения. В численном методе вместо оты- 
скания непрерывных функций времени рассчиты- 
ваются значения искомых функций в дискретные 
моменты времени, а система дифференциальных 
уравнений заменяется тем или иным способом раз- 
ностными уравнениями. В применяемом алгорит- 
ме использована неявная схема Эйлера с пересче- 
том [4] , которая обладает третьим порядком точно- 
сти на шаге и вторым на интервале. 

Полученные системы нелинейных уравнений 
решаются итерационным методом Ньютона [5] на 
каждом временном слое. При этом процесс полу- 
чения системы линейных уравнений, относящихся 
к одной ступени, объединяется с их частичным ре- 
шением. Величины приращений давлений Р Мт , 
Р Тт , Р<) ІШ относящиеся к одной ступени выражают- 
ся через приращения давлений в трассах питания 
ступеней Р І)ш . В результате получаем систему ли- 
нейных уравнений относительно приращений да- 
влений в трассах питания ступеней каскадов. В 
предложенной расчетной модели система линей- 
ных уравнений на каждом шаге итерации решается 
методом исключения Гаусса [4]. 

На каждом временном шаге решение может до- 
полняться расчетом нестационарной гидравлики 
линий МКК, состоящих из подкачивающих ком- 
прессоров, регуляторов давления различных типов 
и других элементов. 

По предлагаемой методике проведены числен- 
ные исследования для различных случаев неста- 
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ционарных возмущений ТС, результаты которых 
позволили сделать следующие выводы: 

• Результаты расчетов, для которых имеются фак- 
тические данные (изменение величин внешних 
потоков, закрытие блоков), позволяют судить о 
достаточной адекватности описанной модели 
поведению реального объекта; 

• Расчетное время распространения гидравличе- 
ских возмущений по ступеням каскадов и ли- 
ниям МКК практически совпадает с имеющи- 
мися опытными данными; 

• Временные характеристики модели позволяют 
оперативно производить расчеты многокаскад- 
ной технологической схемы в режиме реально- 
го времени. 

Заключение 

Предложена математическая модель нестацио- 
нарных гидравлических процессов в промышлен- 
ных центрифужных каскадах, учитывающая отли- 
чительные особенности построения исследуемой 
технологической схемы разделительных произ- 
водств. Математическая модель представляет со- 
бой систему дифференциальных уравнений перво- 
го порядка с заданными начальными условиями и 
позволяет получать решения в режиме реального 
времени при произвольном наборе возмущающих 
воздействий. 

Разработан алгоритм решения уравнений, опи- 
сывающих нестационарные гидравлические про- 
цессы, как в центрифужных каскадах, так и в ли- 
ниях МКК, объединяющих эти каскады. 

Предложенная модель позволяет исследовать не- 
стационарные процессы, возникающие при измене- 
нии величин внешних потоков; отключении внеш- 
них потоков; изменении частоты и отключении пи- 
тающего тока; нарушении вакуумной плотности 
оборудования; изменении точки подачи питания; 
отказах и перенастройке межступенных регуляторов; 
отключении/включении части оборудования; нару- 
шениях в работе оборудования линий МКК; а также 
при различных сочетаниях этих воздействий. 

Модель может быть использована в АСУТС в 
качестве экспертной системы и в компьютерном 
тренажере для подготовки специалистов раздели- 
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